Reaching the hydrodynamic regime in a Bose-Einstein condensate by suppression of 

avalanches 
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We report the realization of a Bose-Einstein condensate (BEC) in the hydrodynamic regime. 
The hydrodynamic regime is reached by evaporative cooling at a relative low density suppressing 
the effect of avalanches. With the suppression of avalanches a BEC containing 120-10® atoms is 
produced. The coUisional opacity can be tuned from the coUisionless regime to a collisional opacity 
of more than 3 by compressing the trap after condensation. In the collisional opaque regime a 
significant heating of the cloud at time scales shorter than half of the radial trap period is measured. 
This is direct proof that the BEC is hydrodynamic. 

PACS numbers: 03.75.Kk,32.80.Pj 



The behavior of excitations in a BEC with energies 
larger than the mean field energy is determined by the 
mean free path of the atoms. Usually the mean free path 
of the atoms is larger than the size of the sample (the col- 
lisionless regime). It would be of great interest to realize 
a BEC in the hydrodynamic regime, where the mean free 
path of the atom is less than the size of the condensate. In 
this situation the properties of the BEC are strongly in- 
fluenced by the inter- atomic collisions. A hydrodynamic 
BEC would give the opportunity to investigate interest- 
ing properties of the condensate, for example, thermal ex- 
citations, heat conduction, shape oscillations, when there 
is only locally thermal equilibrium. 

The transition from the coUisionless to the hydrody- 
namic regime above the BEC transition temperature has 
been studied theoretically Q and experimentally 0] ■ For 
the situation below the BEC transition temperature the- 
oretical discussions are given in Refs. In this 
Letter we will discuss the experimental realization of a 
BEC in the hydrodynamic regime and study the excita- 
tions generated by three-body decay. 

The most obvious way of reaching the hydrodynamic 
regime is creating a large and dense BEC. However, like 
shown in Ref. PI the atom losses will increase strongly 
due to avalanches at such high densities, that are nor- 
mally necessary for entering the hydrodynamic regime. 
This will severely limit the lifetime of the condensate in 
the hydrodynamic regime, as well as the collisional opac- 
ity. Therefore, Schuster et al. Q concluded that the 
collisionally opaque regime can hardly be reached in al- 
kali BEC experiments. In the case of metastable helium 
BEC experiments to reach the hydrodynamic regime has 
so far been unsuccessful H, ■ A second way of realiz- 
ing a hydrodynamic BEC is increasing the cross section 
for the elastic collisions by means of a Feshbach reso- 
nance. The increase in the cross section results in a large 
loss rate which makes it an unsuitable approach for 
achieving a hydrodynamic BEC. Notable exceptions are 
the BEC's of molecules consisting of fermions tuned close 
to the unitarity limit [ij . 

In this Letter we will show that it is possible to enter 



the hydrodynamic regime by following the first approach 
with a strong reduction of the effects of the problems 
mentioned above. This is done by evaporative cooling 
a cloud of atoms to the BEC temperature in a decom- 
pressed trap. The low trap frequencies lead to a rela- 
tive low density of the cloud, which reduces the onset 
of avalanches. The suppression of avalanches results in 
our setup in a BEC containing 120-10^ sodium atoms at 
a density of 2.7-10^* atoms/cm'^. The low density sup- 
presses the losses as a result of avalanches, while the large 
physical size of the BEC results in a collisional opac- 
ity close to the hydrodynamic limit. The hydrodynamic 
regime can subsequently be entered deeply by compress- 
ing the trap after BEC is reached. The lifetime of the 
BEC in the hydrodynamic regime is more than 5 s, which 
is an increase with a factor of 25 compared to earlier 
work 0. 

During the evaporation towards the BEC transition 
the cloud will experience losses due to inelastic collisions 
between the atoms. These losses can be divided into 
three categories: one-, two- and three-body collisions, 
determined by the rate constants G'i,2,3- The one-body 
collisions are the collisions between the trapped atoms 
and the background gas. The pressure of our background 
gas (below 10^^^ mbar, the limit of our detector) results 
in a lifetime of 260 s. These losses are dominant at densi- 
ties below 5.9-10^'^ atoms/cm^ or at temperatures above 
1 ^iK. For reaching a large atom number BEC this is an 
important loss process; however, at the densities reached 
in our BEC it is no longer dominant. The loss rate due to 
two-body collisions is given by G2 • n, with G2 = 6-10^^'' 
cw? /s jljl and n is the density of the sample. This col- 
lision rate never exceeds the one- or three-body losses 
in our experiment and will therefore not be considered 
further in this Letter. 

The dominant loss process at the densities needed for 
the BEC transition in our case is three-body collisions. 
The loss rate for this process is given by G3 • -n? with G3 = 
1.1 •10"'^*' cm^/s 13]. The three-body collisions generate 
high energy atoms that will in the coUisionless regime 
before they rethermalize be removed from the trap by 
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FIG. 1: Number of atoms (a) and the temperature (b) as 
a function of the rf-frequency for evaporative cooling with 
respect to the bottom of the trap. To make the decrease in 
the efficiency more clear two lines are drawn through the data. 



FIG. 2: The number of particles during the evaporation in 
a radial decompressed trap (a) and in an axial decompressed 
trap (b). 



evaporation. However, wlien the sample is in the hy- 
drodynamic regime, the high energy atom will collide 
with other atoms before it leaves the trap resulting in 
an avalanche. In that case a large fraction of the energy 
generated by a three-body decay is released into the BEC 
resulting in a large loss rate of atoms and a strong heat- 
ing of the sample. By contrast, in the coUisionless regime 
only 3 atoms are lost without heating the sample. This 
striking difference in thermal properties makes the hy- 
drodynamic regime very interesting, but extremely hard 
to reach. 

The apparatus is described elsewhere in more detail 
[l3.ll5| . We start the experiment with a magneto-optical 
trap (IVIOT) containing 2.0-10^'' sodium atoms. After 
spin polarizing in a high magnetic field 1.440^" atoms are 
loaded in a cloverleaf magnetic trap (MT) . The magnetic 
field gradient of 118 Gauss/cm in the radial direction and 
the curvature field in the axial direction of 42 G/cm^ lead 
to trap frequencies in the axial and radial direction of 
— 16 Hz and Vp = 99 Hz at a field minimum of 3.4 G. 
We can cool the sample of atoms to degeneracy in 50 
seconds. 

When we measure the number and the temperature of 
the particles during the evaporative cooling, the efficiency 

parameter a. — -^j^ decreases at rf-frequencies below 200 

kHz (Fig.^. The density at this point in the evaporation 
is (7± 1) • 10^'^ atoms/cm^, which is close to the density, 
where the three-body losses become the dominant loss 
process. Despite the large loss rate we can still condense 
35-10^ atoms at a density of 4.0-10^'* atoms/cm^. 

The coUisional opacity k is given by k = (nr) cr^ , with 
(Ts = 87ra~ the s-wave cross section, a the scattering 
length, and (nr) the average column density. Under the 
conditions of Fig.^K =1.3, which means that the mean 
free path of an atom in the condensate is less than its 
size making the BEC hydrodynamic. However, for a thor- 
ough investigation of the hydrodynamics of the BEC it is 
necessary to enter the hydrodynamic regime even deeper. 
Therefore, a BEC containing more atoms is needed. This 



is realized by suppressing the losses due to three-body 
collisions during the evaporation. Since the three-body 
losses scale quadratically with density, this suppression is 
achieved by reducing the density. In order to achieve this, 
the trap is decompressed after 42 seconds of evaporation, 
which corresponds with an rf-frequency of 2.04 MHz with 
respect to the bottom of the trap. The decompression is 
done adiabatically in 2 s. Before this point three-body 
losses are not the limiting loss process because the density 
is always below 2-10^'^ atoms/cm^ yielding a loss rate of 
4-10^'* s^^. For the decompression the trap frequencies 
can be lowered in either axial or radial direction. 

In Fig. [5^ the number of particles during the evap- 
orative cooling is given for a radial decompressed trap. 
The radial trap frequency is in this case lowered from 
99 Hz to 50 Hz resulting in an average trap frequency 
w = ^ujr'-^r'-^z — 34 Hz. Combined with the results of the 
temperature as a function of rf-frequency, we observe no 
decrease in efficiency. The evaporation in this trap results 
in a BEC containing 60T0^ atoms at a density of 2.7T0^^ 
atoms/cm'^. When the decompression takes place in the 
axial direction (16 Hz — > 4 Hz) instead of the radial di- 
rection the efficiency is even higher (Fig.l^b). Note, that 
this decompression results in the same average trap fre- 
quency Lj as above and thus the same decrease in density. 
The evaporation in the axial decompressed trap leads to 
120T0^ atoms in the BEC. This is as far as we know the 
largest condensate starting from an optical trap. Note, 
that we loose less than a factor of 120 in the number 
of atoms during the evaporation indicating an efficient 
suppression of the losses. The increase in the efficiency 
parameter a indicated in Fig. by the dotted lines sug- 
gests that an earlier decompression would even lead to 
more condensed atoms. However, by carefully optimiz- 
ing the sequence it turned out that the decompression 
after 42 s of evaporation leads to the largest number of 
atoms. 

In Fig. |31 the results of a more quantitative analysis 
of decompression are given. In this figure the number 
of atoms in the BEC is plotted for different decompres- 
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FIG. 3: The number of particles in the BEG as a function of 
the axial trap frequency with = 34 Hz. 

sion scenarios. The radial trap frequency is adjusted in 
such way that for every measurement the average trap 
frequency is identical ((D = 34 Hz). The decompression 
discussed above in the axial and radial direction represent 
the two extremes in this figure. The number of particles 
shows a monotone increasing behavior with decreasing 
axial trap frequency. In all cases, the initial tempera- 
ture, density and therefore, three-body collision rate are 
identical. The only difference is the aspect ratio of the 
cloud. In the radial decompressed trap the coUisional 
opacity k is larger, so that atoms undergoing three-body 
collisions are more likely to generate avalanches. Fur- 
thermore, a large n can hamper the evaporative cool- 
ing. Either way, these measurements are a clear sign 
that we have reached the hydrodynamic regime. In the 
axial decompressed trap the avalanches are suppressed, 
until BEC is reached. At the BEC transition the coUi- 
sional opacity has increased to 0.8 resulting in a BEC at 
the cross-over between the collisionless and the hydrody- 
namic regime. Since the suppression of avalanches in this 
case is already effective at 8 Hz, further decompression 
to 4 Hz does not lead to a strong increase in the number 
of particles in the BEC. 

Above we have described how a large atom number 
BEC with K close to 1 is realized by the suppression of 
avalanches. After BEC transition is reached k can be 
tuned from the collisionless regime to a value of 3.2 by 
compressing in 100 ms the trap in axial direction. In 
the remainder of this Letter the adjustability of k is used 
for two kind of experiments. First, the influence of the 
avalanches on the loss rate in investigated. Second, a 
direct proof of the hydrodynamics is presented for k — 
3.2 by studying the heating due to avalanches. 

Previous work 7] has shown that the lifetime of a 
BEC in the hydrodynamic regime is severely limited by 
avalanches. In this experiment the avalanches increases 
the G3 coefficient by a factor of 8. We have measured the 
G3 coefficient as a function of k (Fig. 0)). For each value 
of K the G3 coefficient is determined from the lifetime of 
the BEC. As seen from Fig.01 the G3 coefficient increases 
with more than two orders of magnitude by increasing n. 
At collisional opacities below 0.8 the loss rate is deter- 
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FIG. 4: The increase in G3 coefficient as a function of the 
collisional opacity k. The determined G3 coefficient is divided 
by the G3 coefficient in the collisionless regime fT3| . The solid 
line is a guide to the eye. 



mined by the G3 coefficient in the collisionless regime, 
indicating that there are no avalanches. This is slightly 
above the predicted value of 0.69 for the onset of the 
avalanches [3 ■ This small discrepancy was already men- 
tioned by Streed et al. The difference between the 
predicted value and the measurements indicates that the 
transition from the collisionless regime to the hydrody- 
namic regime is not as sharp as expected. The energy of 
a three-body decay is in the order of 3000 Tc, so it is ex- 
pected that the G3 coefficient will increase even further 
with increasing k. In our experiment we are not able to 
compress the MT further or to condense more atoms, so 
we could not determine the upper limit of G3 . 

When the trap is fully compressed after reaching the 
BEC transition n — 3.2 is reached at a density of 6.6-10^'' 
atoms/cm'^. The increased density due to the compres- 
sion and the effect of the avalanches results in a reduction 
of the lifetime to below 20 ms. The lifetime is sufficient 
for several experiments and can be increased by tuning 
the BEC to a slightly lower collisional opacity. 

In Fig.|31the lifetime of the BEC in the hydrodynamic 
regime is studied in more detail by measuring the num- 
ber of particles in the BEC as a function of the storage 
time for a hydrodynamic BEC with k = 1.7 (circles), and 
a BEC in the collisionless regime with k = 0.8 (squares). 
During the storage time an rf-shield is applied at 100 kHz 
generating a trap depth of 5 /xK to avoid strong heating. 
The lifetime of the condensate in the hydrodynamic and 
collisionless regime is 10 and 20 s, respectively. The hor- 
izontal lines indicate the number of particles needed for 
K = 1.0. Note, that the crossover between the hydro- 
dynamic and collisionless regime is different for the two 
situations, because the trap frequencies are different. We 
can see that the hydrodynamic BEC stays in the hydro- 
dynamic regime for more than 5 s, which is an increase 
with more than a factor of 25 compared to earlier work 
01 • The hydrodynamic BEC was produced with radial 
decompression during the evaporation, proving that we 
can evaporative cool in the hydrodynamic regime. Since 
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FIG. 5: The number of particles in the BEG as a function 
of the storage time for a hydrodynamic BEG with k = 1.7 
(circles), and a BEG in the collisionless regime with k = 0.8 
(squares). The horizontal lines indicates k = 1. 
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FIG. 6: The heating of the sample for a strongly hydrody- 
namic BEG with K — 3.2. The solid line is a guide to the 
eye. 



the losses are larger during the coohng the BEC con- 



tains less atoms. Note, that the decay of the hydrody- 
namic BEC is much faster at small time scales indicating 
avalanche enhanced losses. 

The collisional opacity of 3.2 results in a s-wave colli- 
sion probabihty of 0.96. This means that an atom has 
chance of 96 % to collide with another atom, while travel- 
ing from the center to the edge of the condensate. In this 
situation a significant amount of the energy generated in 
a three-body collision is distributed over the BEC. In the 
collisionless regime the high energy atom will be evapo- 
rated by the rf-shield before the energy is released into 
the BEC. Therefore, in the hydrodynamic regime heat- 
ing can occur on a time scale faster than half the radial 
trap period, which is 5 ms in our case. When we measure 
the temperature of the strongly hydrodynamic sample as 
a function of the storage time we see an increase of the 
temperature due to heating of more than 35% within 4 
ms (Fig.ini). This is a direct proof that the BEC is deeply 
in the hydrodynamic regime. Even when an rf-shield is 
applied this heating is not significantly suppressed. 

In this Letter we have described the experimental real- 
ization of a hydrodynamic BEC. This is achieved by sup- 
pressing the avalanches during the evaporation towards 
the BEC transition resulting in a BEC containing 120-10® 
sodium atoms with a collisional opacity of 0.8. The col- 
lisional opacity is low enough to suppress the avalanches 
during generation of the BEC, while it is high enough to 
yield the possibility to enter the hydrodynamic regime 
deeply by compressing the trap, after the BEC transi- 
tion is reached. The collisional opacity can be tuned from 
the collisionless to the hydrodynamic regime with a colli- 
sional opacity of 3.2 by axial compression. Furthermore, 
the lifetime of the BEC in the hydrodynamic regime is 
more than 5 s, which makes it a very good starting point 
for further research on the hydrodynamic behavior of the 
BEC. 
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